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Abstract: The effects of tlie ion temperature are discussed in a two-ion electron 
plasma and for a model applicable to the oscillating sheath theory that has recently been 
much in the focus of research. The differences between the fluid and kinetic models 
are pointed out, as well as the differences between the approximative kinetic description 
(which involves the expansion of the plasma dispersion function), and the exact kinetic 
description. It is shown that the approximative kinetic description, first, can not describe 
the additional acoustic mode which naturally exists in the plasma with an additional 
ion population with a finite temperature, and, second, it yields an inaccurate Landau 
damping of the bulk ion acoustic mode. The reasons for these two failures of the model 
are described. In addition to this, a fluid model is presented that is capable of capturing 
both of these features that are missing in the approximative kinetic description, i.e., 
two (fast and slow) ion acoustic modes, and the corresponding Landau damping of both 
modes. 



PACS numbers: 52.30.Ex; 52.35.Fp; 52.40.Kh 



I. INTRODUCTION 



In the recent past there has been a lot of activity in the domain of theory and ex- 
periments deahng with the sheath in multi-ion plasmas. An earlier work on this issue^ 
provides us with a general expression for the Bohm criterion, while more recent theoretical^ 
and experimental^ works specify that, in a two-ion electron plasma, the sheath velocity of 
the two ion species have the same value that equals the hybrid ion acoustic (lA) velocity. 
Yet, there seems to be no complete consensus about this in the literature. 

The behavior of the ion acoustic mode in the presence of an additional ion species may 
be drastically different compared to its behavior in a single ion plasma.^^^*^ In particular, 
the Landau damping is increased with the addition of a lighter ion specie, which implies the 
possibility for controlling the mode behavior. Physically, this is equivalent to lowering the 
electron temperature and reducing the phase velocity, so that the damping is increased.^ 
Experimental observations of those phenomena are available in the classic Refs. 8, 9, and 
more recently in Ref. 10. We observe that these effects of the additional ions are just the 
opposite of the effects in the two (hot and cold) electron population case.^^'^^ 

We stress that each additional ion specie introduces a new branch of acoustic oscil- 
lations in the system, complemented by various novel physical effects. However, in the 
context of the sheath theory, this domain of behavior of thermal ions has not yet been 
properly explored. This is the subject of the present work. In the case of a plasma with 
two ion species, one obtains two acoustic modes: one fast and one slow mode. As a result, 
the commonly used terms like the 'system sound velocity' appear to be redundant and 
may become inapplicable in most plasmas. 

II. COMPARISON BETWEEN THE COLD AND HOT IONS CASES 
A. Cold ions 

In order to see how thermal (kinetic) effects influence the oscillating plasma sheath, 
we first solve the dispersion equation for the ion acoustic (lA) mode propagating in a 
three component plasma consisting of electrons, and two cold ion species denoted by a 
and b: 



Here, similar to Ref. 1, we have k = kXde, <^ is normalized to ujpa, vjo = fjo/cso denotes 
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the directed sheath currents/velocities of the two ion species, j = a, b, and c^^ = uTe/nia. 

In the hmit of large wave-lengths and for negligible sheath currents, with the given 
normalization, the dispersion properties of the hybrid bulk lA mode are described by 

uj^ = P[l+ maUbo/ {nibnao)]. (2) 

In physical units it reads u^/k'^ = = {nao/neo)KTe/ma + {nbo/neo)KTe/mb, where Vs 
denotes the hybrid lA sound speed in plasmas with two ion species, called 'the system 
sound speed' in the literature.^ 

The equation ([1]) for cold ions is solved numerically. First, similar to Ref. 1, we set 
Vao = 0.009, Vbo = 0.011, which should correspond to the situation deep inside the plasma, 
far from the sheath, and we also choose rria = 10 m^, and solve for Ubo = 0.2 in terms of 
k. The result (for the ion acoustic mode) is presented in Fig. 1 with a dashed line, showing 
the tendency for the usual frequency saturation for large wave-numbers. Going to larger 
values of k has no sense because of physical reasons. Two additional complex-conjugate 
solutions for the sheath current driven modes, that also follow from Eq. ([1]), are about 2 
orders of magnitude lower and are not of interest here (see Section III) . 

B. Hot ions 

Now, for comparison, keeping the ion thermal effects, we re-derive the dispersion 
equation using the kinetic theory. It reads 

A{uj, k) = l + J2HJk%lJ [1 - Z{uj^o/kv,^)] = 0. (3) 

a 

Here, a = e,a,b, cUao = to — kvao, and Z{x) = [x/(27r)-'^/^] / dy exp^—y"^ /2) / (x — y) is 
the plasma dispersion function, where x = uiao/ {kvj.a), y = v/v^a- For non-streaming 
electrons and in the limit 

^Ta.fe < ^/fc < ^^Te, (4) 

the standard expansions for Z are used, and the general dispersion equation ([3]) in that 
case becomes: 
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\kv^„ ^ 3kv^ _,(^/2)i/2_5^ exp \~ulJ{2k\lJ 



1^ + ^ -'('^/2)V=^exp [-.1,1 (2kV,,)]] = 0, (5) 



In the appropriate limits its real part yields Eq. (jl]), while the wave damping is given ap- 
proximately by 7 ~ —Im/\{uj^k) /[dRe/\{uj^k) / duj\t^=i^^.. Here, ImA{uj, k) and Re A {uj, k) 
denote the imaginary and real parts of Eq. ([5]). 

We have solved Eq. ([5]) numerically by setting Ta = Tb = Te/15, and for the same 
Vjo, Ubo as above. This result is also given in Fig. 1 (by full and dotted lines). The Landau 
damping has the maximum at A; ^ 1 where IjI/uj — 0.25. Note that here and further in 
the text we present the absolute value of the Landau damping. 

Observe a considerable change in the mode frequency. This may be of even greater 
importance than the damping because, compared to the cold ion limit, the phase speed 
is increased, which should be taken into account in modeling the Bohm criterion. 

The presence of the second ion species affects both the mode frequency and the Landau 
damping, as shown in Fig. 2, where the hybrid lA mode frequency (normalized to ujpa) is 
presented in terms of the number density of the species b, and for Ta = = Te/30, and 
= 0.1 (cf. Fig 1). Here, in the beginning the Landau damping is increased by increasing 
the amount of the ion species 6, but this goes only up to rib/ria of about 6 percent when 
\'^\/uj 0.11. Note that for n^o = we have I7I = 2- 10~^ (in units of ujpa), i.e., an 80 times 
lower damping compared with the case when riho/nao = 0.06! Clearly, the presence of the 
additional ion species significantly affects the damping of the mode. This is because the 
wave-particle exchange of energy is easily facilitated when additional ions are available in 
the system. 

For hotter ions (i.e., Ta = = Tg/lO), the ratio 171/^; is increased, with the maximum 
~ 0.22 at ribo/nao = 0.2, as seen from Fig. 3. Here, the Landau damping grows with the 
addition of more 6-ions, yet this increase in the damping saturates for the number density 
Ub reaching about 20 percent. Note that for Ubo = it is I7I = 4 - 10^'^, which is 12.5 lower 
compared with the case at Ubo/nao = 0.2. 

Hence, the second ion species may drastically increase the Landau damping of the 
hybrid lA mode, and we see also that its relative effect is more pronounced for lower ion 
temperatures. 

III. SHEATH-CURRENT DRIVEN MODES 

As is well known, ^'^ the presence of directed fluxes of plasma components implies the 
presence of extra collective modes. In our case these are the ion sheath-current (SC) 
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driven, complex-conjugate solutions with frequencies of the order of ~ kvjo, and with a 
growth rate that is approximately one order of magnitude lower. Their presence is obvious 
from the form of Eq. ([T]) and they are discussed below. 

Increasing the directed ion velocities, the lA mode behavior is changed. This is visible 
from Fig. 4, where the velocities Vjq are multiplied by a factor 100, becoming close to the 
values expected for the Bohm criterion, and for rif, = 0.2na, Ta = Tf, = Tg/lO. Here, the 
normally negative solution for the hybrid lA mode (that is not discussed above in Fig. 1 
as being of no interest) is so strongly Doppler shifted that it becomes positive (line b). 
The absolute value of the Landau damping (the full line in Fig. 5) reaches its maximum 
at ~ 0.7 when |7|/t<;,, = 0.16. 

On the other hand, the growing part of the two complex-conjugate SC solutions, in 
the k domain presented in Fig. 5, has a growth rate that is about one order of magnitude 
below the mode frequency. It can easily be shown that for even larger values of > 15 
the instability vanishes and the real part of the SC mode frequency splits into two modes. 
However, this implies wavelengths below the electron Debye length and the model has no 
sense in that limit. 



IV. ADDITIONAL ACOUSTIC MODE 



Apart from the demonstrated Landau damping and the modification of the frequency, 
the additional hot ion species further imply additional branches of acoustic oscillations. 
This can only be seen by solving the general dispersion equation ([3]) numerically, like 
in Refs. 7, 9. Formally, the additional acoustic branches appear in the limit when the 
terms k kept finite in comparison with u'^. Therefore, in view of the condition 

(jl]), these additional acoustic branches are absent in the previous analysis and can not be 
deduced from Eq. ([5]). Hence, the approximative kinetic approach describes the Landau 
damping but it is unable to describe the additional acoustic branch. However, a recently 
developed fluid theory^"^ easily captures both of these features within the fluid theory. For 
Boltzmannian electrons and with the help of the Poisson equation, the derivation of the 
dispersion equation is straightforward yielding: 

1 + ^ = ^ + . (6) 

Here, Uji = ujo + ifijok"^. In the cold ions limit and for the given normalization, Eq. ([6]) 
becomes identical to Eq. ([1]). The term fijQ = Xvgj/ i^n'^dj), where vl^ = c^j + v^j, follows 



from the ion momentum equation of the form 



{d/dt + VjoV)vji = —{qj/mj)d(f)i/dx — {KTj/njo)dnji/dx + fijod'^vji/dx'^ , 

where it has been introduced to describe the Landau damping. Such a fluid model has first 
been used in Ref. 15, and more recently in Refs. 14, 16. It is convenient because it allows 
the use of the fluid theory, where the expressions can be analyzed more easily, as compared 
to the general kinetic expression ([3]) containing the plasma dispersion function expressed 
through the integral Z{x) = [a;/(27r)^/^] / (i?/exp(— ?/^/2)/(x — y). The term dj = 5j/\ in 
fijo gives the ratio of the Landau attenuation length 6j and wave-length A. It is chosen 
in such a way that it is independent of the wavelength and the plasma density, and it 
depends on the ion temperature in a prescribed way. As is shown in Refs. 14, 16, such 
a fluid model for the intrinsically kinetic Landau damping is, first, much more accurate 
than the approximative kinetic Landau damping obtained after the expansion of Z{x) 
over the parameter x. Second, it gives the same damping as the exact kinetic Landau 
term. Third, it yields the analytical expression for the acoustic modes in two-ion electron 
plasmas with different temperatures of the two ion species [the alternative is numerically 
solving Eq. ([3])]. 

All these features follow after adopting the following expression: ^^'^^ 

dj = Sj/X ^ 0.2750708 + 0.0420737789rj + 0.0890326rJ-0.011785rj' + 0.0012186r/. (7) 

Simple derivations^^ for an electron ion plasma yield the modeled fluid Landau damping 
as: 

|7/|/^ = l/(27rrf). (8) 
The approximative^^ kinetic Landau damping is: 

= {n/Sy/' [{mjm^y/' + r(3 + r)V2 exp [-(3 + r)/2]] . (9) 

Finally, the exact kinetic Landau damping can be obtained only numerically. However, 
using the graph of the exact kinetic Landau damping from Ref. 11 one finds^^ that it may 
be expressed by the following polynomial: 

\-fex\/uJ = 0.681874545 - 0.369763643r + 0.0934588855r2 

- 0.01203427r3 + 0.0007523967677r'* - 0.000018r^ (10) 
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The above given statements about the accuracy of the expressions (18|9|) can directly be 
checked by plotting the expressions flHl fTOj) in terms of r, see Fig. 2 from Ref. 14. 

It appears that, in the case of a two- ion electron plasma, the behavior of two acoustic 
modes is determined by the relation between Csa and Vrb- More details are available in 
Refs. 9, 14. Eq. ([6]) gives two (fast and slow) ion acoustic modes (that appear due to the 
fact that ions have some temperature), together with the corresponding Landau damping. 
As will be demonstrated below, this raises the question about the appropriateness of the 
'system sound velocity', that is commonly used in the literature. 

Using the same normalization as earlier, Eq. ([6]) becomes 

^ + ^ = { ~ ^^aO\AW^) ^ - kVaO^Jrieo/naO + i[kVsa/ (71"rfa)]\AW^ 

-k'^Taneo/iTenao)] 

+ ~~{{'^~ kvbo^Jneo/riao^ UJ - kvbo\/n^o/nao + i[kvsb/ {i^db)]^ rieo/nao 

-k'^nianeoTb/imbnaoTe)^ . (11) 



riao rrib 



Here, v'^j = {cj.j + v'^j)/c1^. For comparison with previous cases, we solve Eq. (ITT]) by 
taking the same parameters as in Fig. 1, i.e., n^o = 0.2nao, Vao = 0.009, Vbo = 0.011, 
Vsa = 1.03, Vsb = 3.26, rria = lOm^. We have also chosen Tg = 11600 K, Ta = Tb = T^/ld, 
and nib = 4mp so that Cga = 1547 m/s. The two acoustic modes are presented in Fig. 6. 
The given sheath currents make negligible frequency shifts in the two acoustic modes. 
The corresponding Landau damping of the two modes, presented in Fig. 7, shows an 
obvious difference as compared to Fig. 1. It is, first, lower by magnitude and, second, 
it shows no decrease for larger values of k. Hence, the decrease and saturation of the 
Landau damping, seen in previous figures, is clearly only due to the expansion of the 
plasma dispersion function, i.e., an artefact of the approximation used, and may not have 
any relevance to the real physical systems. 

Similarly, when Ubo is increased the frequency and the Landau damping of the upper 
(fast) lA mode is increased. However, for the slow lA mode these both parameters are 
reduced when n^o is increased. This is checked by setting k = 0.3 and for other parameters 
as above. The results are presented in Fig. 8 (for the two acoustic frequencies), and Fig. 9 
for the absolute value of the Landau damping (multiplied by 10^). This behavior is in 
agreement with results existing in the literature. ^'^^ 
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Setting the ion temperature to a larger value, i.e., Ta = Ti, = Te/2, the frequencies are 
increased by about a factor 1.5. However, the Landau damping is increased by a factor 
100. This is presented in Figs. 10 and 11. It is seen also that the damping of the slow 
mode is now increased for larger values of Ubo- 

V. SUMMARY 

To summarize, the ion thermal effects are investigated in the sheath theory in multi- 
component plasmas containing two ion species and electrons. Typically, this implies using 
the kinetic theory and we have demonstrated the differences between some results existing 
in the literature (obtained without ion thermal effects), and those obtained in the present 
analysis when the ion temperature effects are taken into account. The most important 
additional effect which follows from the finite ion temperature is the Landau damping, 
that is discussed in detail. However, the standard analysis implies the expansion of the 
plasma dispersion function, and as a result a) some phenomena related to ion thermal 
effects are missed in the procedure, and b) the mode behavior in some limits may show 
features that are not physical and that are only due to the mentioned expansion. The 
most obvious example of a) is the presence of an additional acoustic mode that can be seen 
only by solving the general kinetic dispersion equation ([3]) numerically.'''^ As for b), such 
an artefact of the expansion is the saturation and decrease of the Landau damping, as 
demonstrated in Sees. 2, 3. In the case of hot ions, the expansion of the plasma dispersion 
function is not justified and an alternative approach is needed. We have demonstrated in 
Sec. 4 that such a method exists within the framework of fluid theory. It is simple and 
suitable for analytical work, and at the same time is capable of capturing both effects, the 
additional ion acoustic mode and the Landau damping. In our earlier works^^'^^ we have 
shown that, in fact, quantitatively it describes the Landau damping much more accurately 
than the analysis which follows after the expansion of the plasma dispersion function. 

We conclude that more care is needed whenever the standard expansion of the plasma 
dispersion function is used. This may directly be seen by comparing figures from Sections 2 
and 3 on one side, and those from Sec. 4 on the other side. We stress also that even 
relatively small ion thermal effects imply the presence of an extra (slow) acoustic branch 
of ion oscillations, so that in fact two acoustic modes propagate in the plasma. The 
examples given in Sec. 4 show that the frequencies of the two modes are not very distant 
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even for relatively cool ions (in the given case the ion temperature of only about 1/15 of 
the electron temperature), yet the two modes have rather different behavior. Therefore, 
terms like 'system sound velocity' may become inapplicable for most plasmas. 

A possible extension of the results presented here could be the case of a two-ion 
electron plasma with one negative ion specie, and in particular the pair-ion electron plasma 
produced recently in the laboratory^^ and discussed in the literature very extensively. 
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Figure captions: 

Fig. 1: The positive solutions for the hybrid lA mode frequency (normaUzed to ujpa) for 
hot (full line), and cold ions (dashed line), and the absolute value of the normalized 
Landau damping I7I (dotted line). 

Fig. 2: Full and dashed lines: the positive frequency (normalized to ujpa) of the hybrid 
lA mode in terms of the number density of the second ion species, for k = 0.1. 
Dotted line: the absolute value of the normalized Landau damping |7|. 

Fig. 3: The same as in Fig. 2 but for hotter ions Ta^Tb = Te/10. 

Fig. 4: The positive lA mode (hne a), and its normally negative counterpart, here 
strongly positively Doppler-shifted (hne b). Dashed line: the sheath-current (SC) 
mode. Here k = kXde- 

Fig. 5: Full line: the absolute value of the I A Landau damping I7I (corresponding to 
Fig 4) ; dashed hne: the corresponding growth rate of the SC mode. 

Fig. 6: Two ion acoustic modes in two- ion electron plasma with = = Tg/15 (other 
parameters given in the text). 

Fig. 7: The absolute Landau damping I7I = \'y\/ujpa (multiphed by 10^) of the two lA 
modes from Fig. 6. 

Fig. 8: The frequencies of the two ion acoustic modes in terms of the number density of 
the ion species 6 for = = Tg/lS. 

Fig. 9: The absolute Landau damping I7I = \^\/ujpa (multiphed by 10^) of the two lA 
modes from Fig. 8. 

Fig. 10: The frequencies of the two ion acoustic modes in terms of the number density 
of the ion species b for — — Tg/2. 

Fig. 11: The absolute Landau damping I7I = \'y\/u}pa of the two lA modes from Fig. 10 
(forT„ = Tb = Te/2). 



12 



3 








Landau damping 

O M ^ 



O ■ 



N)|- 




0.2 

bO aO 



0.4 



0.35 



0.30 



0.25 



0.20 



0.15 



0.10 



• fast lA mode 
■ slow lA mode 



0.2 

bO aO 



0.4 



